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Exceptional points are singularities of the spectrum and 
wave functions which occur in connection with level repul- 
sion. They are accessible in experiments using dissipative 
systems. It is shown that the wave function at an exceptional 
point is one specific superposition of two wave functions which 
are themselves specified by the exceptional point. The phase 
relation of this superposition brings about a chirality which 
should be detectable in an experiment. 

PACS numbers: 03.65.Bz, 02.30Dk, 05.45Gg 



Level repulsion is a well known pattern in virtually all 
aspects of quantum mechanics. It states that the levels 
of a selfadjoint Hamiltonian H generically do not cross as 
a function of a parameter A on which H(\) depends [Q. 
Its importance is particularly pronounced in the realm 
of quantum chaos . The connection between excep- 
tional points Q and the occurrence of level repulsion has 
been discussed in || . 

An exceptional point (EP) is a value A c of the pa- 
rameter A, where two of the eigenvalues Ek of H are 
equal to each other - say E U (X C ) = E v+ i{\ c ) - but where 
the space of the corresponding eigenvectors is only one- 
dimensional. We call this a coalescence of the eigenvalues 
and the eigenfunctions \ip u ), It is well known that 

this cannot occur for a selfadjoint Hamiltonian, where 
E v = E u+ \ entails a two-dimensional space of eigenvec- 
tors, in which case the phenomenon is called a degener- 
acy. 

Consider 



oo 

E v+1 (\) = E°„ + Ve s (-v/A-A c ) 



(2) 



H — Hq + XHi , 



(1) 



where Hq,Hi are real and symmetric N x N matrices, 
and let A be a complex number. Then H is a complex 
symmetric matrix. At an EP there is always a singularity 
- namely a branch point - in the spectrum Ek(X) and the 
eigenfunctions \tpk(X)}- The spectrum consists of the val- 
ues that one analytic function assumes on N Riemannian 
sheets. The sheets are connected by N(N—1) square root 
branch points, the EP's. If an EP - connecting E v and 
E u+ i - occurs sufficiently close to the real A-axis, the two 
levels undergo a level repulsion as A sweeps over the real 
axis in the vicinity of the EP. Conversely, when two lev- 
els undergo repulsion, there is always a nearby EP, where 
the expansions 



E y {\) = E° v + J2 e s (V*-*cY 



.1=1 



exist with a finite radius of convergence. 

Three major results have been shown in Q and exper- 
imentally verified in Q when an EP is encircled in the 
complex A-plane: 

1. The two energy levels E v and E„ + i connected at 
the EP are interchanged by a complete turn in the 
A-plane. 

2. The two wave functions \ij)J) and are not 
just interchanged like their eigenenergies but one of 
them undergoes a change of sign. In other words, a 
complete loop in the A-plane leads to {ip v , Vv+i} - * 
{— Vv+ij Vv}- As an immediate consequence we 
conclude: (i) the EP is a fourth order branch point 
for the wave functions and (ii) different directions of 
going through the loop yield different phase behav- 
ior. In fact, encircling the EP a second time in the 
same direction we obtain {—if> u , — Vv+i} while the 
next loop yields {Vv+i, — Vv} an d omv the fourth 
loop restores the original pair {tpu, ipu+i}- It fol- 
lows that the opposite direction yields after the first 
loop what is obtained after three loops in the for- 
mer case. 

3. The behavior of the two energy levels is distinctly 
different when a path in the A-plane is taken below 
or above an EP. In one of the cases, the two levels 
avoid each other while their widths cross, in the 
other case, the two levels cross while their widths 
avoid each other. 

In the topological structure of an EP has been 
shown in the laboratory to be a physical reality. In the 
present paper, we focus attention upon the wave func- 
tion at the EP. We show that the chiral behaviour that 
appears under item (2) above, is an intrinsic property 
of an EP. We argue that this chiral behaviour should be 
detectable in a suitable experiment. 

Recall that for A — > A c , one has |i/v(A)) — ► |"0ep) and 
|VV-H.(A)) — > |"0ep) for the two coalescing wave functions. 
We mention that all the other TV — 2 wave functions are 
regular at a given EP. 

Since H of Eq. (jl|) is not selfadjoint for complex A, the 
right hand eigenvectors \ipk) arc different from the left 
hand eigenvectors (ipk\- Both systems together form a 
biorthogonal basis, i.e. the completeness relation reads 
for A ^ A c 



1 



E 



\^k){^k 



= 1. 



Recall that 



(3) 



(4) 



Due to the symmetric form of H the left hand eigenvec- 
tor {tp\ is just the complex conjugate of its right hand 
partner. Hence, in the Dirac notation, the (complex) 
components of the row vector (ip | coincide with the com- 
ponents of the column vector From Eq.(Q) it follows 
that 



EP 







(5) 



since the orthogonality holds identically in A and thus 
in particular at A — A c , when j — v and j' = v + 1. 
As a consequence, the inverse of the biorthogonal norm 
(V'felV'fe) that appears in Eq.(|^), does not exist at A = A c 
for k = v, v + 1 . 

For a two-dimensional space, N = 2, one concludes 
from Eq.(||) that IV'ep) has the form 



EP 



±1 
1 



(6) 



Nowhere in the foregoing, a basis has been laid down 
with respect to which the coefficients of the vector are to 
be taken. In fact, (^) remains true under all orthogonal 
transformations of a given basis, even complex orthog- 
onal ones. These are the transformations that conserve 
the symmetry of H which we consider. Hence, in every 
basis with respect to which H is symmetric, I^ep) will 
have the form (|^). In particular, there is no orthogo- 

i \ I —i 



nal transformation that maps 



onto 



1 



Every 



\ipEp) is therefore either 



For illustrative purpose we consider a two level model 
in detail. We stress, however, and below explicitly elabo- 
rate that even an infinite dimensional problem is, in the 
vicinity of an EP, locally equivalent to a two dimensional 
problem. 

Consider 



H 



ei 
e 2 



with 



U(<f>) 



xu [ ^ ) U T 

uj 2 



cos cp — sin ( 
sin d> cos 6 



(7) 



(8) 



where the angle <f> and the energies ek,^k, k — 1,2 are 
real. The eigenvalues are 

El AX) = ei + e2 + t {uJl + UJ2) ±R, 0) 



where 



R= < (^) 2 + (^r^) 2 



(10) 



+ 2^( £ i - £ 2)(^i ~ U2) cos 20 



1/2 



The two levels coalesce when i?(A) vanishes. This hap- 
pens at 



A 



£l ~ £2 
UJl - UJ 2 



exp(±2i<; 



(11) 



Note that for zero coupling {cj> = 0) the two branch points 
cancel each other and a genuine degeneracy occurs with 
the well known properties of a diabolic point 0. 

The eigenfunctions can be parametrised by the com- 
plex angle 9 as follows 



Here, 9 is related to the parameters in Eq.(0) via 

tan 6»(A) = A(wi - u 2 ) sin 2cj>/ (13) 
{Ei(\) - E 2 {\)+e 1 - e 2 + A(cji - w 2 )cos20). 

The eigenvectors of Eq.([l2]) are normalised in the 
biorthogonal sense 



= ! for fc = l,2, A^A C 



(14) 



An explicit calculation shows that the coefficients of the 
wave functions Eq.(|lJ) diverge at the EP. Inserting A^ 
into Eq.(|13|) we obtain 



tan 6^ 



This implies for A — > Xf 
cos 9 ± 



sm t 



00 



(15) 



While the completeness relation Eq.(||) is obeyed for 
those values of A which do not coincide with an EP, the 
set of eigenfunctions is incomplete at an EP. At the EP, 
the eigenfunctions (O) coalesce for A — * A^ as 



|^i (A)> F 1 



1^2 (A)) - F 2 



±i 



±i 
1 



(16) 



Here, the factors F\ t 2 depend on A; in fact they diverge 
for A — > A c . 

So far we have established that, in the two dimensional 
model, the wave function at the EP has, up to a com- 
plex factor, a strictly prescribed form: the ratio of the 
components is +i at A+ and —i at A~. This holds in 



2 



any basis and irrespective of the parameters ei, u>i and 
(j). Before we turn to the physical relevance of this result 
we discuss a higher dimensional situation as this would 
usually prevail in an experimental set-up. 

In higher dimensions the wave function at the EP will 
of course no longer have the simple form of Eq.(|l6|), in 
fact |"0ep) has then N components. But using the com- 
pleteness relation of Eq. (0) we expand 



with 



where 



and 



|lM=5> fc (A)| Xfc (A)) 
fc 

Cfc = (%*#ep), 



(17) 



IXfc-; 



(Xk 



(ipk\4>k) 



(fAfcl^fc 



which ensures biorthogonal normalisation. If this expan- 
sion is used in close vicinity of an EP, where \ip v ) an d 
are about to coalesce, it is obvious that only the 
terms in Eq. ( |l7| ) with fc — v and fc = v + 1 make substan- 
tial contributions. In fact, all Ck vanish when A — > A c as 
follows from the orthogonality for fc =/= v, v + 1 and from 
Eq.(||) for fc = v, v + 1. However, the vanishing numera- 
tors for fc = v, v + 1 are compensated by the vanishing de- 
nominators with the result that in the limit A — ► A c only 
the terms with k = v and fc = v + 1 survive. This result 
implies that the ./V-dimensional vector IV'ep) is basically 
a superposition of only the two (iV-dimensional) vectors 
|Vv(A)) and |i/v+i(A)); the closer A is to A c the more cor- 
rect is the statement. In other words, with regards to the 
EP, the iV-dimensional problem can be locally simulated 
by a two-dimensional problem. From Eq.(|l6|) we thus 
conclude that 

= +i or = — i (18) 



must hold in the vicinity of A c but independent of A 
within this vicinity. 

A more explicit analytic consideration shows how this 
astounding and important result comes about. We de- 
note the components of \4>ep) by {xk}, fc = 1, . . . , N and 
Ik x k = 0. If A is near to A c the 



recall (Eq.(§)) that £ 
components of the unnormalised |Vv(A)) can be chosen as 
{x k + d k } with d k = a fe VA - A c + 0(A-A c ), fc = 1, . . . , N 
and some constants dk being of no interest here. The 
components of the unnormalised (^„+i(A)| must there- 
fore, to lowest order in the dk, have the form {xk — dk}. 
To lowest order in the dk we obtain 



(i>u\ifep, 
(ipv+ilifep, 



= 2y^x k d k 

fc 

fc 

= ^ x kdk 

k 

= — / J Xkdk 



(19) 
(20) 
(21) 
(22) 



from which the statement of Eq.(|18|) immediately follows. 

The local reduction - in the vicinity of an EP - of the 
full iV-dimensional problem to an effective two dimen- 
sional problem is now achieved by the two-dimensional 
matrix h = /iq + \h\ with the matrix elements 



(Mii' = (Xj\H \Xj>) 



J, J 



= v,v + l 



(23) 
(24) 



using the relevant state vectors \Xu) and |x^+i)- In Fig.l 
we display two different but typical examples to demon- 
strate how efficiently the procedure works. The eigen- 
values of hg and hi yield the effective values of the tj 
and u)j as used in Eq.(^). The effective coupling an- 
gle 4> is obtained from the eigenvectors of hi in the ba- 
sis where ho is diagonal (note that ho from Eq.(|23|) is 
not a priori diagonal). The straight lines in Fig.l are 
the lines tj + XtOj which correspond to the effective un- 
perturbed lines (</> = 0); switching on (f> to the calcu- 
lated value yields an almost exact approximation of the 
iV-dimensional problem by the effective two-dimensional 
problem. For each level repulsion, i.e. for each EP, the 
procedure has to be carried out from the outset. The 
example of Fig.l is based on a random ten dimensional 
case. 
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FIG. 1. Two-dimensional approximations of the EP asso- 
ciated with level repulsions. The top drawing displays a sec- 
tion of a ten-dimensional problem. The drawings in the mid- 
dle and at the bottom are blow-ups of the encircled areas of 
the top. The straight lines are explained in the text. The dis- 
tinction between the exact and the effective two-dimensional 
problem is within the line thickness for the curved lines in the 
middle and bottom drawing. 

To summarise: an EP is locally equivalent to a two- 
dimensional problem. Knowing all parameters of the ef- 
fective two-dimensional problem we know, from Eqs.jl^) 
and (|l8|), the specific superposition of the wave function 
at the EP in terms of those wave functions which coalesce 
at the EP. We find the relation 



|-0ep) = +i\Xv) 
for A+ 

IV'ep) = ~i\Xu) 
for A7 



lx»+i) 

e v ~ £ ^+l +2i<t> 
UJ V — U) u +l 

\Xu+l) 

— £u+l c -2i<j> 



(25) 



(26) 



In a higher dimensional problem the quantities tj , tOj and 
4> arc effective quantities as defined above. 

In an experimental situation like a microwave res- 
onator, the phase factor +i means that the time depen- 
dent wave function \x v ) has a leading phase of a quarter 
of a full period with respect to |Xi/+i). For the phase — i 
the wave is lagging by the same amount. This should 



be detectable ||. In the particular case, where the two 
wave functions can be associated with two independent 
linear polarisations, the wave function at the EP would 
then be an elliptic or circular wave with a definite chiral- 
ity. We note that a similar observation has been made 
in 0] for the treatment of damped acoustic waves in a 
solid medium. If the two wave functions can be associ- 
ated with different parities, the superposition again has 
a definite chirality. 

We conclude that a definite chirality is associated with 
each EP. In a high dimensional problem one expects a 
random occurrence of a particular chiral behaviour just 
as the random occurrence of the associated level repul- 
sions. Note that, in an experiment, only those EP are 
accessible which have a negative imaginary part of the 
eigenenergy. Depending on the effective values of the 
€j , ojj and 4>, these points may lie in the upper or lower 
A-planc. 
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